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We present a minimal physical model for the flagellar motor that enables bacteria to swim. Our 
model explains the experimentally measured torque-speed relationship of the proton-driven E. coli 
motor at various pH and temperature conditions. In particular, the dramatic drop of torque at high 
rotation speeds (the "knee") is shown to arise from saturation of the proton flux. Moreover, we 
show that shot noise in the proton current dominates the diffusion of motor rotation at low loads. 
This suggests a new way to probe the discreteness of the energy source, analogous to measurements 
of charge quantization in superconducting tunnel junctions. 



The bacterial flagellar motor is a molecular machine 
that rotates a helical filament and thereby powers the 
swimming of bacteria like E. coli [T^]. Motor rotation is 
typically driven by H + ions that generate torque by pass- 
ing into the cell via the motor, down an electro-chemical 
gradient called the proton motive force (PMF) . Although 
much work has been devoted to understanding proton 
translocation and its coupling to torque generation, bio- 
chemical details are lacking and many questions remain 
unanswered. An important one is whether ion translo- 
cation is cooperative, i.e. whether protons translocate 
individually or in groups. Here, we present a minimal 
physical model for torque generation (Fig. [TJ that not 
only explains a variety of previous experimental obser- 
vations, but also suggests a way to measure the cooper- 
ativity of proton translocation. Specifically, the model 
predicts that at low loads, motor diffusion is dominated 
by proton shot noise with a strong (quadratic) depen- 
dence on proton cooperativity. 

The flagellar motor operates with near-perfect effi- 
ciency at low speeds [TU]. As the speed is increased, 
e.g. by reducing the load, the torque and efficiency ini- 
tially remain high — the "plateau" of the torque-speed re- 
lationship (TSR) — and then drop abruptly at a "knee" 
(c/. Fig. [2]). This knee occurs at higher speeds as 
temperature is increased. Despite much experimental 
011 G3J [[6] and modeling progress U GU HS1 ED] , 
the origin of the knee is still poorly understood. In [T]l2"0] 
the cause of the knee was argued to be the gating of pro- 
ton translocation by the relative position between stator 
and rotor. In [3], a detailed model of motor kinetics was 
proposed to explain the observation that motor speed is 
independent of the number of stators at low loads [21] . 
In our model, proton translocation, which is assumed to 
be the rate-limiting step, is modeled by a barrier crossing 
event. The knee in the TSR then arises from the kinet- 
ically limited rate of proton translocation. Importantly, 
our model fully incorporates proton thermodynamics and 
yields the separate dependence of the TSR on the elec- 
trical and chemical parts of the PMF. 

Three ingredients underlie our model: (i) Each torque- 
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FIG. 1: Schematic model of the bacterial flagellar motor. 
Left: The passage of a proton, or possibly a group of protons, 
through a torque generating unit (a MotA/B stator — only one 
stator of about 10 is shown) causes a protein spring to stretch 
to its next attachment site, represented by circles, on the ro- 
tor. Right: To translocate a proton must pass through an 
external gate, over a barrier, and finally through an internal 
gate, with all the steps assumed to be reversible. The net 
energy difference driving proton translocation is the electrical 
potential energy, —eAip, minus the work, t88, necessary to 
stretch the protein spring by 89, where r is the torque applied 
by the spring to the rotor. 



generating unit (MotA/B stator) contributes indepen- 
dently and additively to the total torque, in agreement 
with experimental observations (Fig. ^p) [21 . (ii) The 
torque from each MotA/B stator is applied to the rotor 
by a protein spring. Proton translocation into the cell 
causes the stretching of a protein spring to its next at- 
tachment site (Fig.[l]) [HUH]. This assumption enforces 
the tight coupling between proton current and rotation 
speed [10] , (Hi) Assuming a cooperativity index of n, 
translocation occurs through three reversible steps: first 
n protons load into an external gate, then all n cross an 
energy barrier to an internal gate, and finally all n are 
released into the cell. The barrier crossing event is the 
rate-limiting step. The external and internal gates are 
necessary to explain the non-linear dependence of the 
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TSR on proton concentrations (Fig. |2p,). (We define the 
distance between two attachment sites as n89, so that 
the average displacement per proton is 86.) 

The external and internal gates are assumed to be 
in fast equilibrium with the external and internal pro- 
ton concentrations, respectively. Their dissociation con- 
stants are denoted by K cxt and Ki nt , so that the oc- 



cupancy of the external gate is H™ xt /(K r c 



where iJ cx t is the proton concentration outside the cell, 
and similarly for the internal gate. The energy dif- 
ference between the internal and external gates is n x 
[eAtp + t89 + k B Tlog(K in t/K cxt )}. Aip is the trans- 
membrane electric potential (i/>int — i/'cxt) and t 89 is 
the work necessary to stretch the protein spring. When 
Aip — and t = 0, the energy barrier per pro- 
ton to inward translocation is UP„, and the barrier 



to outward translocation is U® nt , 



with V? - C/ O ut = 



kBTlog(Ki nt /K cx t)- In general, some fraction a of the 
electric potential and some fraction /3 of the work con- 
tribute to the inward barrier, so that the barriers per 
proton to inward and outward translocations are, re- 
spectively, Ui n = U® n + aeAip + fir 66 and U out = 
Uout — (1 — a)eAip — (1 — 0)t86. For simplicity we have 
neglected any dependence of -Kint,ext on Aip and r. Then 
the rate of inward proton translocations is: 
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Jo is a kinetic constant (in Hz), and the other prefactors 
represent the occupancies of the external and internal 
gates. The outward rate J ou t is given by a similar ex- 
pression, so that the net inward proton flux is: 
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Ap := Aip + (k B T/e)log(H int /H cxt ) is the PMF com- 
posed of the electrical and chemical potential differences. 
Its value is approximately —150 mV in normal condi- 
tions. To account for the data, we assume that the 
height of the barrier may depend on temperature, and 

we expand the prefactor to linear order in temperature: 
J e- u °ni T V^T = j e»(T-To) ; where Tq = 17 7 o Cj from 

Fig. [2j>, is chosen as a reference temperature. 

Rotation is then described by coupled stochastic equa- 
tions for the angular position of the rotor 6 and the 
stretching of the protein springs i — 1, . . . , N, where N 
is the number of stators, each exerting a torque n on the 
rotor: 
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v is the frictional drag coefficient of the load, and fc(r) is 
the spring constant of the (possibly non-Hookean) pro- 
tein springs. The spring constant need not be speci- 
fied as none of the observables computed below depend 
on it. In the second equation, each spring relaxes as 
the rotor moves (back-reaction of the rotor onto stators, 
—d6/dt), but gets restreched by proton translocations 
(Cl(n) '■= (Jin — Jout)S9). £(t) is the thermal noise on 
the load and satisfies the Einstein relation: (£(£)£(<')) = 
2(fegT '/ 'i/)S(t—t'). £J is the shot noise at each stator i, due 
to the randomness of proton translocation events. To ob- 
tain the average speed u> and torque per stator r = (rj), 
we solve Eqs. (3j4l at steady state {dn/dt — 0) and in 
the absence of noise, yielding: uo — d6/dt — Nt/v and 
n = t, with fl(r) = lo. Along with Eq. (|2]), this gives a 
closed set of equations from which we obtain the TSR. 
Note that the resulting expression depends separately on 
the electrical and chemical potential for protons. At stall 
(y — * oo) the system is in equilibrium. The energy nec- 
essary for a protein spring to move to its next attach- 
ment site, which is proportional to the torque r it exerts 
on the rotor, is matched by the PMF, t59 + eAp = 0. 
Consequently the total torque grows linearly with the 
PMF, Nt = —NeAp/59, and the efficiency near stall is 
w 100%, in agreement with experiments [5} IT2]. 

Our model with no cooperativity (n = 1) can fit all 
existing measured TSR of the E. coli motor. Some our 
model's parameters are fixed properties of the motor and 
thus are fit by single values: (86 = 4.6°, K int = 1.2 • 10~ 8 , 
JC ex t = 2 • 10~ 7 , J = 670 Hz, a = 0.2, (3 = 0.078, rj = 
0.11 K _1 ) while others depend on conditions (T, i?j n t ext) 
Aip, number of stators N), and may or may not have 
been measured in the experiments. 

Fig. [2^ shows fits of TSRs measured under various pH 
conditions 15J. The electric potential Aip was not mea- 
sured and so was used as a fitting parameter for each set 
of pH conditions. Our fit indicates that \Aip\ increases 
with pH (Fig. |3|, consistently with previous measure- 
ments m ng. 

Electrorotation experiments [5] have been used to 
apply an external torque on the load via an oscillating 
field. When the motor is driven backwards (upper-left 
quadrant of the TSR), the internal torque is approxi- 
mately equal to its stall value up to speeds of — 100 Hz [3] . 
When the motor is driven to speeds larger then the max- 
imum operating speed (lower-right quadrant of theTSR) , 
the motor resists rotation, resulting in a negative internal 
torque. In this regime, the slope of the TSR remains ap- 
proximately the same as for positive torques beyond the 
knee [2]. Our model agrees with measurements in both 
regimes. The absence of a barrier to backward rotation 
follows from the reversibility of proton translocation. For 
negative torques, the model predicts an inflection of the 
TSR, as seen for the red curve in Fig. [2^,. 

Within our model, the proton flux is limited by the 
loading of the external and internal gates and by the 
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FIG. 2: The torque-speed relationship (TSR) of the E. colt 
flagellar motor, a. Rotation speed is measured with beads of 
various loads (data points) attached to a flagellar stub, un- 
der different pH conditions (colors) [15) : solid curves: model 
fits. b. Rotation speed measured with 0.25-0.36 /im diameter 
beads attached to flagellar stubs for a wide range of viscosi- 
ties, at three different temperatures [J; circles: experimental 
data, solid curves: model fits. Inset: Total torque vs. speed 
from [16] . normalized by the number of stators. Data collapse 
indicates that stators contribute independently and additively 
to the total torque. Solid curve: model TSR for a single sta- 
tor using the same parameters as in the main figure, with the 
temperature fit as 21°C and the stall torque fit as 300 pN.nm. 



barrier crossing. This limitation on flux accounts for the 
knee of the TSR. In Fig. [2^i, the position of the knee 
strongly depends on pH values: as the internal proton 
concentration increases, the internal gate gets saturated, 
preventing protons from translocating inwards, and thus 
limiting the stretching of stator springs and the applied 
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FIG. 3: Electric potential difference versus external pH, 
as obtained from the fit to data from [K] (Fig. [2^i) using 
Aip = Ap — (fesT/e) log(//int/// cxt ), and assuming that the 
maximum torque (2500 pN.nm) is reached at Ap = —150 
mV. The apparent linear dependence of Atp on external pH 
is consistent with previous data [8] [13] . 



torque. 

Our model also accounts quantitatively for measure- 
ments of the TSR at different temperatures (Fig.[2]3) [4 , 
as well as for measurements with different numbers of 
stators (Inset of Fig. |2]d) [16]. As a test, we compared 
the model to earlier measurements at different tempera- 
tures [5] (data reported in [4]), with no additional fitting 
parameters, and found excellent agreement (Fig. [4}. 

Experiments by Gabel and Berg (Fig. [5]) [5] have been 
interpreted to imply that the rotation speed is propor- 
tional to PMF, even at high speeds beyond the knee of 
the TSR. Our model predicts that speed is proportional 
to PMF at low speeds, in the plateau region of the TSR. 
However at high speeds, the torque is limited by the pro- 
ton flux, and therefore both torque and speed grow sub- 
linearly with PMF. Nevertheless, our model (with n = 1) 
is fully consistent with the measurements reported in [B] . 
Experimentally, speeds were simultaneously recorded for 
two motors of the same cell, one rotating the cell itself 
(high load) and the other rotating a small polystyrene 
bead (low load), as shown schematically in the lower in- 
set of Fig. [5| These two loads correspond to the two 
dashed lines in the upper inset of Fig. [5] As cells were 
de-energized by the introduction of a respiratory poison, 
the PMF Ap regressed from —150 mV to 0, and the mo- 
tors slowed, with the two speeds approximately propor- 
tional to each other, even at low temperature where the 
low-load, high-speed motor was in the kinetically-limited 
regime. We fitted the data for each cell, using N/v ce \\ 
and N/v bead as free parameters, and assuming that dur- 
ing de-energization the electric and chemical parts of the 
PMF regressed in fixed ratio to each other. The model 
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FIG. 5: Relationship between the high-speed (low-load) and 
low-speed (high-load) regimes of the E. colt flagellar motor, 
as the PMF varies from —150 mV to 0, at temperature T = 
16.2°C. Symbols: Experimental data from individual cells [BJ. 
Solid curves: model fits with the same parameters as in Fig. [2] 
Dotted curve: the model predicts loss of proportionality at 
very low loads. Lower inset: schematic of the experiment [6j. 
Upper inset: model TSR for PMFs of -150, -100, and -50 
mV; the dashed lines show the low and high load lines of the 
cell represented by in the main panel, and the dotted line 
is the load line corresponding to the dotted curve in the main 
panel. 



FIG. 4: The torque-speed relationship at various tempera- 
tures. Comparison between data from [2] and model pre- 
diction (solid curves) with the same parameter values as in 
Fig. [2Jd (no additional fitting parameters). Note that the tem- 
perature of the top panel (11.2°C) is below the range used to 
fit parameters (15.8-22.7°C). 



fits are consistent with the data in the considered pa- 
rameter regime, both at 24°C (Fig. [sj and at 16.2°C 
(Fig. |5j. However, a further reduction of the low load 
would be predicted to lead to a strong deviation from 
proportionality (dotted line in Fig. [5j . 

Adding cooperativity (n > 1) to the model still cap- 
tures the general shape of the TSR, as well as its tem- 
perature dependence (Fig. [6]), but agrees poorly with the 
observed pH dependence (Fig.[7|, as it implies a stronger 
dependence of proton current on ion concentrations. The 
model with n > 1 is consistent with Gabel and Berg's 
data for most cells (Fig. |9| , although it breaks down for 
the two cells with the largest high speeds. 

Although n = 1 appears to best explain the data, 
measurements of average rotation speeds do not allow 
us to discriminate with certainty the proton cooperativ- 
ity. In contrast, our model predicts that diffusion of the 
rotor angle at low loads should depend very strongly 
on proton cooperativity. The discrete nature of pro- 
ton translocations implies the existence of shot noise 



£J{t) in the stretching of the protein spring at each sta- 
tor i. We can approximate the shot noise as Gaus- 
sian white noise: (£f (i)ff(i')) = 2D shot S(t - if), with 
Ashot = (1/2) (Jin + J ut)nS9 2 (this approximation is 
valid for times much larger than Jr 1 and Jout)- Solving 
Eqs. (3j4l including shot noise £J{t) and thermal noise 
£(i), we find an exact expression for the effective diffu- 
sion coefficient of the rotor angle: 
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where N is the number of stators and /i(r) 
— (cifi/dr)" 1 is minus the local slope of the TSR. Fig 
shows the effective diffusion coefficient as a function of 
motor speed for different proton cooperativities. Param- 
eters were chosen so that the speed at zero torque is 200 
Hz. 

At high loads (y ^> /.i), diffusion is entirely due to ther- 
mal noise: D e g kBT/v. However, at low loads (v <C 
n), diffusion is dominated by shot noise: D e ff « D s ^ ot /N. 
In fact, the thermal noise is completely suppressed in the 
low-load limit: e.g., a small thermally induced backward 
jump in rotor angle causes the stretching of all springs, 
which then rapidly pull the rotor forward, thus canceling 
the jump. Notice that in the low-load limit, the shot- 
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FIG. 6: The torque-speed relationship with proton coopera- 
tivity n = 1, ... ,4. Left panels: model fit to data from [4]. 
Right panels: comparison between model prediction and elec- 
trorotation data from [2: with the same parameter values as 
in the left panels. 
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FIG. 7: Model fits of the torque-speed relationships from |15] 
with proton cooperativity n — 2. 



noise contribution is inversely proportional to the num- 
ber of stators. Intuitively, a small jump in the angular 
stretch A#,; of one protein spring ultimately only causes 
the rotor to move A9i/N because the rotor is equally cou- 
pled to all N stators. The variance per jump is therefore 
(A8i/N) 2 and with TV independent stators, the resulting 
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FIG. 8: Relationship between the high speed (low drag) and 
low speed (high drag) regimes of the bacterial flagellar motor, 
as the proton motive force varies from —150 to mV, at tem- 
perature T — 24° C. The experimental data from individual 
cells [6] are represented by symbols. Our fits with no proton 
cooperativity (n = 1) are the solid curves. 
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FIG. 9: The same high speed versus low speed data as Fig. [5] 
but model fits (solid curves) obtained with proton coopera- 
tivity n = 2. 



diffusion scales as 1/N. 

Eq. ^ could be used to infer n experimentally. In 
the low-load limit, we have J m J ou t, and therefore 
n s=s 2D e g/(u)56). 89 can be determined by measuring the 
torque per stator and the PMF at stall: 86 = —eAp/r. 

Our analysis of rotor diffusion suggests a novel exper- 
imental test to investigate the cooperativity of proton 
translocation. Some rotational diffusion measurements 
have already been made [El [18], but not in the regime of 
very low load, where shot noise is expected to dominate. 
Although we have derived the expression for diffusion in 
the specific framework of our minimal spring model, the 
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FIG. 10: Effective diffusion as a function of rotation speed 
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same approach is generalizable to more detailed models 
of the bacterial flagellar motor. 

We thank Avigdor Eldar and Anirvan Sengupta for 
helpful suggestions. T.M. was supported by a Human 
Frontier Science Program fellowship, and N.S.W. by Na- 
tional Institutes of Health Grant R01 GM082938. 



[1] Bai, F., C.-J. Lo, R. M. Berry, and J. Xing, 2009, Biophys 
J 96(8), 3154. 

[2] Berg, H. C, and L. Turner, 1993, Biophys J 65(5), 2201. 
[3] Berry, R. M., and H. C. Berg, 1999, Biophys J 76(1 Pt 
1), 580. 

[4] Chen, X., and H. C. Berg, 2000, Biophys J 78(2), 1036. 
[5] Fung, D. C, and H. C. Berg, 1995, Nature 375(6534), 
809. 

[6] Gabel, C. V., and H. C. Berg, 2003, Proc Natl Acad Sci 

USA 100(15), 8748. 
[7] Iwazawa, J., Y. Imae, and S. Kobayasi, 1993, Biophys J 

64(3), 925. 

[8] Lo, C.-J., M. C. Leake, T. Pilizota, and R. M. Berry, 

2007, Biophys J 93(1), 294. 
[9] Meacci, G., and Y. Tu, 2009, Proc Natl Acad Sci USA 

106(10), 3746. 

[10] Meister, M., and H. C. Berg, 1987, Biophys J 52(3), 413. 
[11] Meister, M., S. R. Caplan, and H. C. Berg, 1989, Biophys 
J 55(5), 905. 

[12] Meister, M., G. Lowe, and H. C. Berg, 1987, Cell 49(5), 
643. 

[13] Minamino, T., Y. Imae, F. Oosawa, Y. Kobayashi, and 

K. Oosawa, 2003, J Bacteriol 185(4), 1190. 
[14] Mora, T., H. Yu, Y. Sowa, and N. S. Wingreen, 2009, 

arXiv physics. bio-ph. 
[15] Nakamura, S., N. Kami-ike, J. ichi P Yokota, S. Kudo, 

T. Minamino, and K. Namba, 2009, J Mol Biol 386(2), 

332. 




